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Abstract 

Many extensions of the standard model predict heavy metastable 
particles which may be modeled as solitons (skyrmions of the Higgs 
field), relating their particle number to a winding number. Previous 
work has shown that the electroweak interactions admit processes in 
which these solitons decay, violating standard model baryon number. 
We motivate the hypothesis that baryon-number-violating decay is a 
generic outcome of collisions between these heavy particles. We do so 
by exploring a 2+1 dimensional theory which also possesses metastable 
skyrmions. We use relaxation techniques to determine the size, shape 
and energy of static solitons in their ground state. These solitons 
could decay by quantum mechanical tunneling. Classically, they are 
metastable: only a finite excitation energy is required to induce their 
decay. We attempt to induce soliton decay in a classical simulation 
by colliding pairs of solitons. We analyze the collision of solitons 
with varying inherent stabilities and varying incident velocities and 
orientations. Our results suggest that winding- number violating decay 
is a generic outcome of collisions. All that is required is sufficient 
(not necessarily very large) incident velocity; no fine-tuning of initial 
conditions is required. 



1 Introduction 



1.1 Motivation 

Many extensions to the standard model which involve strong dynamics at 
the electroweak scale include new heavy particles which have been modeled 
as solitons. The simplest model within which such particles can be analyzed 
is the standard electroweak theory with the Higgs boson mass mn taken to 
infinity and with a Skyrme term added to the Higgs sector. With these 
modifications, the Higgs sector supports a classically stable soliton whose 
mass is of order the weak scale, typically a few TeV. [0] 

To understand how solitons arise, note that in the absence of the 
weak gauge interactions, the Higgs sector of the standard model is a 
four-component scalar field theory in which a global 0(4) symmetry is 
spontaneously broken to 0(3), with vacuum manifold S^. In the tuh — ^ oo 
limit, the dynamics is that of an 0(4) nonlinear sigma model. Field 
configurations are maps from three dimensional space onto S^, and the 
solitons (skyrmions) are configurations which carry the associated winding 
number. The winding number is topological and soliton number is conserved. 

Gauging the weak interactions changes the picture qualitatively because 
the winding number of the Higgs field is not invariant under large gauge 
transformations. This means that a soliton can either be described as 
a skyrmion of the Higgs field with gauge field = or, equivalently, 
as a topologically trivial Higgs field configuration with a suitably chosen 
nonvanishing A^. The latter description makes manifest the fact that there 
are sequences of gauge and Higgs field configurations, beginning with a soliton 
and ending with a vacuum configuration, such that all configurations in the 
sequence have finite energy. This means that the soliton is only metastable: 
it is separated from the vacuum only by a finite energy barrier and can 
decay quantum mechanically by tunneling. [|], ^, |^, |^ Or, the soliton can be 
kicked over the barrier if it is supplied with energy. The process in which an 
electroweak soliton is hit with a classical gauge field pulse (a coherent state of 
ly-bosons) and caused to decay has been analyzed numerically. It is even 
possible to find a limiting case of the theory in which the quantum mechanical 
cross-section for a process in which a soliton is struck by a single VT-boson 
and induced to decay can be calculated analytically. [0] In any process in 
which a soliton is destroyed, one net baryon and one net lepton from each 
standard model generation is anomalously produced. 0] 
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Electroweak solitons have also been studied in the electroweak theory 
with finite Higgs mass, in which the Higgs sector is a hnear sigma model. @] 
If a Skyrme term is added to the theory, metastable electroweak solitons 
exist if rriH is sufficiently large. In the linear sigma model, the Higgs field 
can vanish at a point in space with only finite cost in energy. The Higgs 
winding number is therefore not topological even in the absence of gauge 
interactions. This means that in a world with gauge interactions and a finite 
Higgs mass, there are two ways for solitons to decay: either via nontrivial 
gauge field dynamics, as sketched in the previous paragraph, or via the Higgs 
field itself simply unwinding.]^ 

The metastable electroweak soliton is an intriguing object to study. And 
yet, it is not found in the standard electroweak theory where the Higgs sector 
is a linear sigma model with no higher derivative terms. The Higgs sector of 
the standard model is best thought of as an effective field theory describing 
the low energy (weak scale) dynamics of the light degrees of freedom in 
some higher energy theory. The simplest examples of higher energy theories 
which feature particles which can be described as electroweak solitons in 
the low energy theory are technicolor theories, in which the technibaryons 
play this role. Regardless of whether the underlying theory is specifically 
a technicolor model, it will introduce all higher derivative terms allowed 
by symmetries, including the Skyrme term, into the Lagrangian of the low 
energy effective theory. If the Higgs boson is discovered to be light (say, with 
mass tjih ^ v = 250 GeV), the correct low energy effective field theory will 
almost certainly not support solitons, regardless of the physics of the higher 
derivative terms. If the Higgs boson is discovered to be heavy, there will be 
some class of appropriate high energy theories whose low energy effective field 
theories, although more complicated than that obtained simply by adding a 
Skyrme term to the standard model, feature metastable electroweak solitons. 
Discovery of the corresponding TeV scale particles would confirm that nature 
chooses such a theory. 

Processes in which two metastable electroweak solitons collide have to 
date not been studied. Our purpose in this paper is to use the analysis of a 
two dimensional toy model which shares some (but not all) of the features 
outlined above to motivate the hypothesis that the generic outcome of such 
collisions may be the destruction of one or both solitons. This suggests (but 
certainly does not demonstrate) that baryon number violation is the generic 
outcome of collisions between two of the TeV scale particles which can be 
modeled as solitons. 
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As a sideline, we note that our numerical methods work equally well 
for describing soliton-soliton and soliton-antisoliton collisions. Our focus 
is on soliton decay in soliton-soliton collisions; we note, however, that the 
numerical simulation of soliton-antisoliton annihilation in the Skyrme model 
is well-known as a difficult numerical problem, plagued with instabilities.^ 
We are able to follow soliton-antisoliton annihilation without difficulty (with 
energy conserved at the part in 10^ level). This suggests that our numerical 
methods — in particular the use of the linear sigma model — may be of 
broad utility when generalized to 3+1 dimensions. 

1.2 Metastable Baby Skyrmions 

Let us now introduce the 2+1 dimensional model whose metastable solitons 
we analyze. The Lagrangian density, which describes the dynamics of a three 
component scalar field = {(j)^, 0^, 0^), is 



^dj-d'^f - '^{dj X dp$) ■ (9> X SV) 
— fi^{v — n ■ (f) — X (j) ■ (f — v^^ 



Here, is a unit vector which we choose to be (0, 0, 1). 

To understand the features of this Lagrangian, it is worth beginning by 
setting = and taking the limit A ^ oo. When //^ = 0, the theory has 
an 0(3) symmetry. For A — > oo, one removes the fourth term from ( |1 . 1|) 
and instead imposes the constraint that ■ = i;^ at all points in space and 
time. Because the field is constrained to take values on a two-sphere of 
radius v, field configurations with fixed boundary conditions at infinity can 
be classified by their winding number 

Q = / eafe0 ■ (5a0 X db(j))d^x = I ■ ((9^0 X dy(f) dxdy , (1.2) 

which is integer-valued and topological: configurations with different winding 
number cannot be continuously deformed into one another. This suggests the 



^See Ref. [Q for classical simulations of skyrmion-skyrmion scattering in the 3+1 
dimensional Skyrme model which report instabilities in the simulation of skyrmion- 
antiskyrmion annihilation; see Ref. ll^ for a discussion of the origin of the instabilities 



and Refs.|ll, b_2| for efforts to overcome them. 
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possibility of soliton solutions to the classical equations of motion. Solitons 



in 2+1-dimensional 0(3) sigma models were first discussed in Ref. [|1^], and 
their quantum field theoretic properties were analyzed in Refs. pTSf . Such 
solitons are often called baby skyrmions [0 because of their similarity to 
3+1-dimensional skyrmions. Although our motivation is the analogy to 3+1- 
dimensional electroweak solitons, we note that baby skyrmions themselves do 
arise in certain 2+1-dimensional electron systems which exhibit the quantum 
hall effect [|T6|, although the Lagrangian used in their description differs from 
that in Eq. (|lTTl) . 

The four-derivative term in the Lagrangian ( |1 . 1| ) is the analogue of the 
Skyrme term. It stabilizes putative solitons against shrinking to arbitrarily 
small size. If we were working in three spatial dimensions, the two-derivative 
term would stabilize putative solitons against growing to arbitrarily large 
size. In two spatial dimensions, however, the two- derivative term cannot play 
this role because its contribution to the energy of a configuration is scale 
invariant. We must therefore introduce a zero- derivative term in order to 
stabilize solitons against growing without bound. Such a term must explicitly 
break the 0(3) symmetry, and therefore has no analogue in 3-|-l-dimensional 
electroweak physics, in which no explicit 0(4) symmetry breaking terms 
are allowed. The particular form of the /i^ term in (|l.lj ) therefore has no 
electroweak motivation; it is analogous to a pion mass term in the 3+1 
dimensional Skyrme model, but this is not relevant to us. This model (with 
/i^ nonzero and A — >■ oo) was considered in Ref. |]T7|, and its solitons have 



been analyzed in detail in Refs. [18, 19|. Similar models, differing only in 



the choice of the explicit symmetry breaking term in the Lagrangian, have 
also been analyzed. pO| 



The soliton mass and size in the theory with Lagrangian (|1.1|) with A = oo 
are given by|T^ 



M,„, = 19.47F 



"^^/o:3T6+"^ 



i?3oi~(3-4)/s:^/^, (1.3) 



with ai and 02 dimensionless constants (independent of k/i) satisfying 
Oi + 02 = 1. The parametric dependence of these results can be understood 
by noting that the energy of a configuration of size R receives contributions of 
order Fk,'^R~^ and F/i^R^ from the first three terms in the Lagrangian 
(|1.1| ) and that, as described above, a soliton is stabilized by the balance 



between the four-derivative k"^ term and the zero- derivative /i^ term. 
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If we stopped here, with A infinite, our sohtons would be absolutely stable, 
rather than metastable. Soliton-soliton collisions have been simulated in this 
theory, but of course the solitons never decay. Once A is finite, the fields 
are allowed to deviate from 0-0 = f^, and the soliton configuration with 
■ = f ^ found previously in the A oo theory may unwind and decay. 
Indeed, we will see that soliton solutions do not exist for A less than some Ac. 
If A > Ac, metastable solitons exist: these solitons are classically stable if left 
unperturbed, but can be induced to decay if supplied with sufficient energy. 
Our goal is to determine whether the means by which the energy is delivered 
is important or whether soliton decay is the result of generic soliton-soliton 
collisions, without finely tuned initial conditions. 

For our purposes, A is the most important parameter in the theory because 
by choosing its value, we control the energy required to make the soliton 
decay and indeed control whether solitons exist in the first place. We are 
not interested in the dependence on the other parameters, and indeed most 
of them can be scaled away. We first set f = 1 by rescaling 0. Next, the 
constant F has units of energy and we henceforth measure energy in units 
such that F = 1. Next, k has units of length and we henceforth measure 
length in units such that k = 1. Note that this means that ^ 7^ 1 in our units, 
but this will not concern us as we only discuss the classical physics of this 
model. We have set the speed of light c = 1 throughout. The parameters 
and A~2 are also length scales in the Lagrangian, and the theory is therefore 
fully specified by the two dimensionless parameters Ak^ = A and ^ik = /i. 
Although results do depend on /i, we are not interested in this dependence, 
and we choose to follow Ref. ITS and set /i^ = 0.1 throughout. Once we 



have chosen units with F = k, = 1 and have chosen to set /i^ = 0.1, then 
Msoi = 19.47 and -Rgoi ~ (3 — 4) in the theory with A = 00. 

In Section 2, we find metastable soliton configurations for finite values of 
A with A > Ac ~ 7.6. We shall see that for all values of A for which solitons 
exist, the soliton mass and size change little from their values at A — *• cx3. 
Although we do not explore their dependence on /i, we expect it would be 
similar to that in ( |1.3| ). In Section 3, we present our results on soliton-soliton 
collisions. We find that soliton decay occurs for incident velocities greater 
than some critical value Vc- We explore how this critical velocity depends 
on A and on the initial impact parameter and relative orientation of the two 
solitons. We find that Vc is less than or of order half the speed of light 
regardless of the relative orientation as long as A < 2Ac and b is less than 



5 



or of order the soliton size. Thus, inducing sohton decay does not require 
specially chosen initial conditions; it is a generic outcome of soliton-soliton 



collisions. We make concluding remarks in Section 4. [21 



It perhaps goes without saying that our model is at best a crude toy model 
for the electroweak physics which motivates our analysis. First, we work in 
2+1 dimensions. Second, in order for the theory to have soliton solutions 
we are forced to include a zero-derivative explicit symmetry breaking term 
not present in the electroweak theory. Third, we do not introduce a gauge 
field. Hence, our solitons can only decay via unwinding the scalar field; 
in the electroweak theory, gauge field dynamics introduces a second decay 
mechanism which has no analogue in our theory. Related to this, our solitons 
are absolutely stable for A = oo, whereas electroweak solitons are metastable 
even for mn = oo. This is perhaps the biggest qualitative difference between 
our model and electroweak physics. Fourth, one may worry that even if an 
analysis along the lines of ours were done in the 3+1 dimensional electroweak 
theory itself, the momenta required would make it impossible to analyze 
soliton decay within the effective theory. This concern may be evaded 
for solitons which are almost unstable: in this circumstance, for example, 
VT-soliton collisions can result in soliton destruction even if the VT-boson 
momentum is small enough that the calculation is controlled. Soliton- 
soliton scattering in our model is far from being a complete analogue of 
the scattering of TeV scale particles which can be modeled as metastable 
electroweak solitons; we nevertheless hope that our central result, namely 
that metastable baby skyrmions in 2+1 dimensions are destroyed in collisions 
with generic initial conditions, motivates future work on baryon number 
violating scattering in this sector. 



2 Finding Static Solitons 

Before we can study soliton-soliton collisions, we must find the metastable 
soliton configurations for different values of A. We do this by looking for 
configurations which minimize the static Hamiltonian i^static at a given A. 
The static Hamiltonian is given by 

-f^static = ~ J d^xCstatic ; (2-1) 

where ^static is the Lagrangian density of ( |1.1| ) with all terms containing time 
derivatives set to zero. 
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We discretize -^static on a square lattice of 125 x 125 points, with the 
spatial separation between points given by Ax = 0.2 (in our units in which 
K = 1). We discretize the two derivative term in the standard fashion, writing 
it as a sum over terms like 

(2.2) 



</)*(a;, y) — 0*(x — Ax, y) 



Ax 

The Skyrme term is trickier to handle, because it involves terms like 

d^(f)^ dy(j)^ 9^0^ dy(j)'^ . 



(2.3) 



We discretize this contribution to the Hamiltonian as a sum over terms like 
' (f)^{x + Ax, y) — (f)^{x — Ax, y)\ [ (j^ (x, y + Ax) — 0^ (x, y — Ax) 



2Ax J V 2Ax 

0^(x + Ax, y) — (t?{x — Ax, y)\ ( 4>^{x, y + Ax) — 0^(x, y — Ax) 



2Ax 



2Ax 



(2.4) 



In this way, we ensure that within each term in the sum over lattice sites, 
all spatial derivatives are centered at the same point in space. Discretizing 
the Hamiltonian in this fashion ensures that discretization errors are of order 
(Ax)2. 

In order to find a soliton, we begin with a guess (which we describe 
momentarily) for the configuration 0(x, y) and perform a numerical 
minimization of the static Hamiltonian using the conjugate gradient method 
of Ref. [Q. (It is important to use a method such as this one, which 
minimizes a function of N variables using computer memory of order 
rather than of order iV^ since we have an iV = 3 x 125 x 125 dimensional 
configuration space.) In order to minimize the energy, the conjugate gradient 
routine needs expressions for the gradient of the energy at any point in our 

dimensional configuration space, with respect to each direction in this 
configuration space. We obtain these expressions by varying the discretized 
-f^static with respect to the 0* at each lattice site. (These expressions will of 
course also appear as the terms with no time derivatives in the dynamical 
equations of motion of Section 3.) 

For A 



oo, soliton solutions can be written in the form 17, IS 



sin /(r) cos 6* 
{r,9) = I sin /(r) sin 
cos /(r) 



(2.5) 
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where /(r) satisfies the following conditions: 

/(O) = TT, (2.6) 
lim /(r) = . (2.7) 

r— >oo 

(We define polar coordinates such that the soliton is centered at r = 0, 6* = 
is the positive y-axis, and 9 increases in a clockwise direction.) Note that 
because of the /i^ term in the Lagrangian which breaks the 0(3) symmetry, 
must point in the (jr" direction at large r. The 0(2) symmetry associated 
with rotations in the ((/)^,0^) plane is not broken in the Lagrangian; in the 
solution, these rotations are mapped onto rotation in the (x, y) plane about 
the soliton center. This configuration is thus a two-dimensional analogue of 
what in three dimensions is called a hedgehog configuration. 

In our search for solitons at finite A, we therefore begin by choosing 
a reasonably large A, namely A = 15, and making an initial guess of the 
form (|2.5|) with f(r) = 7rexp(— r/2). We then run the conjugate gradient 
relaxation algorithm repeatedly, until the change in the energy between 
successive relaxation steps is smaller than one part in 1O^°.0 The soliton 
configuration we find is a hedgehog configuration, as at A — oo. However, 
when A is finite, (f) ■ (f) ^ 1. The soliton we find can be written in the form 

(sin/(r) cos6' \ 
sin /(r) sin ^ (2.8) 
cos fir) ) 

with fir) satisfying the same boundary conditions as above.0 We depict the 
soliton configuration in Fig. 

^As a check, we then used this configuration as an initial condition for the full time- 
dependent dynamical equations of motion described in the next section. The total kinetic 
energy during the time evolution was never more than one part in 10^ of the soliton energy. 
This confirms that the relaxation algorithm has indeed converged to a static solution to 
the full equations of motion. 

■^Note that we could have rewritten the static Hamiltonian in terms of cr(r) and /(r), 
discretized that Hamiltonian in r, and then used a conjugate gradient algorithm to find 
these two functions of r. This would have been less computationally intensive than finding 
(j3^{x, y) as we did. However, the expressions we obtain by varying our static Hamiltonian 
relative to the fields at each lattice site, and indeed the results we obtain for at each 
lattice site in a soliton configuration, are precisely what we need in the next section when 
we analyze soliton-soliton collisions, which are of course not circularly symmetric and so 
cannot be written in terms of (j{r) and ,f{r). 
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Figure 1: /(r), a{r) and the energy density for the sohtons with A = 15 
(sohd curves) and A = 7.7 (dashed curves). 



After obtaining a baby skyrmion at A = 15, wc used the resulting 
configuration as the initial condition for relaxation at A = 14, and so found 
the soliton configuration at this A. We repeated this process step-by-step in A, 
finding solitons for values of A down to A = 8. At A = 7, energy minimization 
led to a configuration with zero energy, instead of to a soliton. We then used 
the A = 8 soliton as an initial configuration for relaxation at A = 7.9, and so 
on down to A = 7.6 where again no soliton was found. We therefore know 
that a stable soliton exists at A = 7.7. It is a logical possibihty that there is 
a stable soliton at A = 7.6 even though our relaxation algorithm did not find 
one. We think this is unlikely, because the soliton configurations which we 
have found at A = 7.7 and A = 7.8 are very similar, and we therefore believe 
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1 9 1 50S 


1 3 
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12 


19.0751 


11 


19.0250 


10 


18.9618 


9 


18.8791 


8 


18.7619 


7.9 


18.7473 


7.8 


18.7309 


7.7 


18.7131 


7.6 


no soliton 



Table 1: Energy of Static Solitons at various Lambdas. 



that the A = 7.7 sohton is a very good starting configuration from which to 
find the A = 7.6 sohton if it existed. We therefore conclude that classically 
stable solitons exist only for A > Ac, with 7.6 < Ac < 7.7. 

In Table I, we give the energies of the solitons which we have found for 
various values of A. In Fig. |l|, we depict the field configuration and energy 
density for the solitons we have obtained for A = 15 and A = 7.7. We note 
that even though A = 7.7 is only just above Ac, the soliton configuration does 
not look very different from that at much larger values of A, and the soliton 
energy is also little changed. Note that the deviation from air) = \ is only at 
most 20% for a soliton with A = 7.7 which is on the edge of instability. The 
central energy density does increase by almost a factor of two as A is reduced 
from 15 to 7.7. Note, however, that the total energy is almost unchanged, 
and actually decreases slightly. The soliton radius decreases as A is reduced 
towards Ac, but does not decrease dramatically. The definition of R^o\ is of 
course somewhat arbitrary; if we take it to be the radius inside which 90% 
of the total energy of the soliton is found, we find -Rsoi = 3.31 for A = 15 and 
i?soi = 2.83 for A = 7.7. 

Although the energy density and \J 4> ■ <p = <J are circularly symmetric, 
the fields 0^ and 0^ in a soliton configuration are not circularly symmetric. If 
we only observed a single static soliton, this would be of no consequence: in a 
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hedgehog configuration, the different possible choices for cfy" and 0^ are related 
simply by rotations in space. However, when we describe a configuration of 
two well-separated solitons in the next Section, the relative angle a between 
their orientations does matter. That is, specifying such a configuration 
requires giving the relative position and velocity of the centers of the two 
solitons and the angle a. The first soliton in such a configuration can be 
mapped onto the second by a translation followed by a rotation by an angle 
a about the soliton center. 



3 Colliding Solitons 

With solitons in hand, we are ready to study what happens when they 
collide. For this purpose, we need discretized equations of motion and 
a numerical algorithm to evolve an initial configuration, now specified by 
0* and (ff at each lattice site, forward in time. We begin by writing a 
discretized Lagrangian which is a function of 0* and 0* at each of the 
lattice sites, at a single time t. We discretize the time-independent terms 
as described in the previous Section. There are no spatial derivatives of 
0* in the Lagrangian, so discretizing terms involving 0' is trivial. We then 
use the Euler-Lagrange procedure on this Lagrangian written in terms of 
3 X 125 X 125 0's and 3 x 125 x 125 0's, and obtain equations of motion 
which specify the 3 x 125 x 125 0's. These equations of motion take the form 
of three coupled linear equations for 0^, 0^ and 0'^ at a given lattice site, 
which are easily solved. We now have an expression for (j)^{t,x,y) written in 
terms of the values of 0* and 0* at lattice sites within two spatial links of 
the site of interest, all at the same time t.0 We are now ready to take a step 
forward in time. 

We evolve the system forward in time using the Runge-Kutta-Feldberg 
algorithm and the adaptive algorithm of Ref. [E^ for choosing the size of 



the time step At. That is, we first use the fifth-order Runge-Kutta-Feldberg 
algorithm to obtain 0* and 0* at time t+At. This fifth-order method is special 
because a rearrangement of the fifth-order function evaluation terms results 
in a fourth-order Runge-Kutta expression. We then have two different 

^Note that because of the way we discretize spatial derivatives in the Lagrangian, 
expressions in the equations of motion with mixed time-space derivatives such as dtdx(f>^ 
end up discretized as [(p^{x + Ax, y) — (jf{x — Ax, y)]/2/S.x. 

■^This hidden fourth-order expression is referred to as an embedded Runge-Kutta formula 
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estimates (fourth order and fifth order) for 0' at t + At at each lattice site, 
and can evaluate the discrepancy between the two estimates for each of the 
3 X 125 X 125 0's and 0's. If the largest discrepancy is larger than a specified 
tolerance, we reject the step and begin anew with a smaller At. We use 
the largest discrepancy to estimate how much At should be reduced. If all 
discrepancies are smaller than the specified tolerance, we accept the result 
of the fifth-order calculation for 0* and 0* at time t + At. After a successful 
step forward in time, we use the largest discrepancy (which must have been 
less than the tolerance since the step forward was accepted) to estimate by 
how much we can safely increase At when we take our next step forward in 
time. In the simulations of collisions which we describe below, the tolerance 
is such that the timestep selected by the adaptive algorithm is approximately 
0.01 < At < 0.05. Note that we do not use conservation of energy as our 
criterion for acceptance or rejection of a step forward in time. This makes it 
fair to use a check of the conservation of energy as an independent measure of 
the accuracy of our evolution algorithm. We do this at various points below. 

We choose fixed boundary conditions, with fixed to its vacuum value 
(0, 0, cTvac) at the boundaries of our 125 x 125 grid, where o"vac solves 

2 2 

("^vac ~ l)cvac = "^vac — 1 + ^ ^r large A. Since the solitons 

have radii of order Rsoi ~ 3, we choose initial conditions with two solitons 
whose centers are a distance 10 apart. We initialize (p by adding these two 
soliton configurations. (That is, we take ^vacuum + (0first soiiton - ^vacuum) + 

— * — * 

(^second soliton - ^vacuum)-) The resulting configuration is not precisely a 
minimum of the static Hamiltonian, but the two solitons are far enough 
apart that this is not a big concern. To obtain a soliton moving with an 
initial speed v in the positive x direction, we simply initialize 

^\x, y) = -v[ct>\x, y) - ct>\x - Ax, y)]/Ax (3.1) 

at time zero. For simplicity, we are using a Galilean boost. This is 
appropriate for w -C 1. When we use this prescription with a velocity at 
which relativistic corrections are becoming important, the initial condition 
we have specified is not the correct Lorentz-boosted, Lorentz-contracted 
soliton. In this circumstance, as the system is evolved forward in time, 
the soliton radiates some energy and quickly settles down to become a 
(correct) relativistic soliton moving with a velocity somewhat less than v. 

due to the fact that it can be obtained with no additional function evaluations. 
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Figure 2: Sequence of snapshots of the energy density during a coUision 

between two sohtons which resuhs in the destruction of both. The grey scale 
indicates energy density. In this simulation, A = 10, the initial velocity of 
each soliton is f = 0.5, the impact parameter is 6 = 0, and the solitons have 
a relative orientation angle a = in the initial configuration. The images are 
at times t — 0,4:, 8, 12, 16, 20. In this and in all subsequent figures showing 
soliton-soliton collisions, each panel shows a 25 x 25 square (in our units in 
which K, = 1) and the initial separation between solitons is 10. The lattice 
spacing is Ax = 0.2. 

For example, when we set v = 0.8 in our Galilean boost prescription for the 
initial condition, we in fact end up with a soliton moving at a speed of 0.61. 

We begin by analyzing collisions between two solitons in the theory with 
A = 10. We choose initial conditions in which both solitons are moving 
(towards each other) with velocity v = 0.25, with zero impact parameter. 
We choose an initial relative orientation angle a = 0, meaning that one 
soliton is obtained from the other by translation without rotation. Previous 
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Figure 3: Left panel: Kinetic, potential and total energies during the soliton- 
soliton collision shown in Fig. with A = 10 and v = 0.5. The topmost curve 
(constant to better than two parts in 10^) is the total energy. Of the other 
two curves, the one that begins low is the kinetic energy, the one that begins 
high is the potential energy, including spatial gradient energy. Right panel: 
Same, during the time-reversed evolution. We reverse the sign of all 0* in 
the final configuration of Fig. 0, and then watch the evolution algorithm 
recreate the initial configuration of Fig. |^. 



work shows that two static solitons with this relative orientation repel each 
other. ||19[ This is consistent with what we find: for low velocities, as for 
example for v = 0.25, the two solitons bounce off each other and return 
whence they came. We now increase v to 0.5. This time, the outcome, 
depicted in Fig. ||, is that the solitons are destroyed in the collision. The 
final state is a cloud of debris, namely small amplitude oscillations of the 
field spreading outwards from the scene of the collision. In Fig. |^, we 
show the kinetic energy, potential energy and total energy for the collision 
shown in Fig. ^. (By "potential energy" we mean the contribution to the 
energy from all those terms in the Hamiltonian with no time derivatives. 
Most of this energy is due to spatial gradients of the fields.) First, we see 
that the total energy is conserved, in fact to better than two parts in 10^. 
The kinetic energy is not zero initially, because the solitons are moving. As 
the solitons approach each other, the kinetic energy decreases. This confirms 
that the interaction is repulsive: the solitons slow down and deform as they 
approach. As the solitons approach each other more closely, at some point 
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their deformation becomes sufficient that they are no longer stable, and they 
fall apart. The resulting outgoing waves have approximately equal kinetic 
and potential energy, as expected for traveling waves. It is quite clear from 
Fig. ^ if it was not already clear from Fig. ^ that the solitons have been 
destroyed. 

As a stringent check of the accuracy of our time evolution algorithm, we 
take the final configuration from our simulation, reverse the sign of 0*, and 
evolve it for the same period of time as we did initially. The second panel 
of Fig. ^ shows the behavior of the energies during this "backwards-in-time" 
evolution. It is clear that the debris reconstitutes itself into two solitons! 
The sequence of snapshots of the energy density looks almost exactly like 
those in Fig. 0, but in the opposite order in time. The discrepancies between 
the energy density in the initial configuration and that in the configuration 
obtained after soliton collision and destruction followed by time-reversed 
evolution and soliton recreation differ by at most 1/40 of the energy density 
at the center of the soliton. The total energy is conserved to better than one 
part in 10^. 

As a further check of the stability of our algorithm, we have also simulated 
soliton-antisoliton annihilation. We obtain an antisoliton configuration from 
a soliton configuration by making the transformation cfp' —0^, equivalent 
to taking 9 — > —9 in ( p^.5D or ( ^.81 ). This turns a hedgehog configuration into 
an anti-hedgehog configuration, and hence yields an antisoliton. We find 
that analyzing soliton-antisoliton collisions using our evolution algorithm is 
no more difficult than analyzing soliton-soliton collisions. We were able to 
follow the annihilation process with energy conserved to better than one part 
in lO'*. Now, with confidence in the accuracy and stability of our evolution 
algorithm, we proceed to analyze the outcome of soliton-soliton collisions 
with a variety of initial conditions. 

We ffist explore how the outcome of a collision depends on A and w, 
keeping the impact parameter 6 = and the relative orientation angle a = 
as above. The results of many simulations are summarized in Fig. ^ We 
discover that for any A, there is a critical velocity Vc below which the solitons 
rebound without decaying, and above which one or both (usually both) 
solitons are destroyed. This critical velocity goes to zero as A — > Ac. As 
A is increased, Vc increases, reaching about half the speed of light for A about 
twice Ac 

We now return to A = 10, f = 0.5, still keeping a = and ask how the 
outcome of a collision depends on the impact parameter h. For h = 2.0, both 
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Figure 4: Outcome of soliton-soliton collisions with different initial velocities 
and different values of the parameter A. All collisions have impact parameter 
6 = and relative orientation angle a = 0. Note that v is the velocity 
parameter in ( |3.1| ). If v is large enough that relativistic effects are significant, 
the actual velocity of the soliton is somewhat less than v. For example, 
V = 0.8 yields a soliton with velocity 0.61. 

solitons decayed into traveling waves, as we found for 6 = above. We show 
the outcome of this collision in Fig. ^. Note that b = 2.0 is a substantial 
impact parameter, comparable to the soliton radius Rgoi — 3. We find that 
the solitons still decay if 6 = 3.2. An impact parameter b = 4.0, however, 
yields a collision which is sufficiently peripheral that the solitons emerge 
intact, deflected from their initial directions of motion by about 45°. We 
can describe our results by saying that the critical velocity Vc above which 
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Figure 5: Snapshots of energy density during a collision between two solitons 
with impact parameter b = 2.0 in the theory with A = 10. The relative 
orientation angle is a = 0. The initial velocity v = 0.5 is large enough that 
the solitons are destroyed. The time between images is 4.0. 

soliton decay is the outcome of the collision increases with increasing impact 
parameter. For 6 = 0, Fig. |^ shows that 0.27 < Vc < 0.3. We now see 
that Vc = 0.5 for a nonzero impact parameter in the range 3.2 < b < 4.0. 
We have also done several more simulations with b = 2.0 and various initial 
velocities, and find that for b = 2.0, the critical velocity is 0.3 < Vc < 0.4. We 
conclude that soliton decay does not require collisions with small or finely- 
tuned impact parameters. Although increasing b from zero increases the 
critical velocity Vc required to destroy the solitons somewhat, it remains easy 
to destroy solitons as long as the impact parameter is less than or comparable 
to the soliton radius. 

All the collisions we have described to this point have had the same 
relative orientation. For a = 0, low velocity collisions yield a rebound, in 
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Figure 6: Snapshots of energy density during a collision between solitons with 
relative orientation a = 180°, impact parameter 6 = 0, and initial velocity 
V = 0.25 in the theory with A = 10. The solitons are not destroyed and 
(eventually) form a classically stable bound state. The time interval between 
images varies: the images are at times t — 0,4:, 8, 12, 16, 20, 24, 28, 34, 42, 50. 
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which each sohton reverses direction, while higher velocity collisions lead to 
soliton destruction. We now consider a collision (with A = 10, f = 0.25 and 
6 = 0) between two solitons with a relative orientation angle a = 180°. That 
is, the second soliton in the initial configuration is obtainable from the first 
by a translation and a 180° rotation. The interaction between static solitons 
with this orientation is known to be attractive.! 19 1 We show the outcome of 



a low velocity collision in Fig. y. The work of Ref. |T8[ reveals that in the 
A oo theory, there is a stable, ring-shaped, soliton with winding number 
2. It appears that the final state of the collision in Fig. ^ will be a soliton 
of this form, although it will differ in its details from that of Ref. [0 since 
A is finite. What we observe in Fig. ^ is that the incident solitons at first 
scatter by 90°, but then do not escape to infinity. They fall back upon one 
another, and rescatter by 90°. There are small outgoing ripples at late time, 
but they have too little energy density to be visible in Fig. ^. We expect 
that were we to run the simulation for a long time, in a big enough box that 
outgoing ripples never return, we would see repeated 90° scatterings, with 
the solitons escaping less and less far away each time, all the while radiating 
small outgoing ripples, and eventually settling down to become the static, 
ring-shaped configuration. 

As we increase the incident velocity, we find that for v > Vc with 
0.43 < Vc < 0.48, the outcome of the collision is soliton destruction rather 
than 90° degree scattering followed by the formation of a bound state. 
We show an example of collision induced decay in a collision with relative 
orientation a = 180° in Fig. ^. Note that the critical velocity above 
which soliton destruction is the outcome is somewhat larger than, but still 
comparable to, that we found previously for a = 0. We have not mapped 
out Vc vs. A for the a = 180° orientation as we did in Fig. ^, but we expect 
that the figure would be qualitatively similar. One new feature, though, 
would be that at large A there would be two different outcomes possible for 
collisions with v < Vc- bound state formation (for low enough v) and 90° 
scattering followed by the escape of the two intact solitons to infinity (for 
larger v which is still less than f^). At A = 10, we do not find any velocities 
for which 90° scattering followed by escape occurs. It must occur at larger A, 
since it certainly occurs at large enough velocities for A ^ cxo, when fc — 1- 

The collision shown in Fig. |^ is an example of a simulation in which the 
initial velocity {v = 0.5 in this case) is only just above the critical velocity 
(0.43 < Vc < 0.48 in this case). In this circumstance, what we generically 
observe is that the solitons scatter, separate a little, but are sufficiently 
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Figure 7: Snapshots of energy density during a collision between two solitons 
with relative orientation a = 180° in the theory with A = 10. The impact 
parameter is 6 = 0. The initial velocity is large enough {v = 0.5) that the 
two solitons decay. The time between each image is 4.0. 



distorted as a result of the scattering that after separating a little they fall 
apart. We observe this phenomenon also at a = 0, except in this case the 
solitons scatter by bouncing back in the direction whence they came, then 
separate a little, and then fall apart. At velocities which are somewhat larger 
than Vc, as for example in the collision shown in Fig. 0, we find that soliton 
destruction occurs more promptly, during the initial collision. 

We now consider collisions between solitons with a relative orientation 
angle a = 90°, still with A = 10 and 6 = 0. For this relative orientation, there 
is no force between static solitons. |jl9| We find the same possible outcomes 
as we did for a = 180°. As a function of increasing velocity, the outcome of 
a collision is either capture to form the ring-shaped bound state, or soliton 
destruction. (Again, scattering by an angle of 90° followed by the escape of 
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two intact solitons would be a possibility at larger A.) The critical velocity 
above which soliton decay occurs is 0.25 < Vc < 0.3. 

4 Concluding Remarks 

We have analyzed soliton-soliton collisions in a 2 + 1-dimensional theory 
with metastable baby skyrmion solutions. We find classically stable soliton 
solutions for values of the parameter A which are larger than Ac ~ 7.6. These 
solitons are prevented from decaying by a finite energy barrier and so can 
decay if supplied with sufficient energy, for example in a collision with a 
second soliton. We have mapped out the space of initial conditions under 
which the outcome of a soliton-soliton collision is the destruction of one 
or both solitons. We find that soliton decay results whenever two solitons 
collide with an incident velocity greater than some Vc- This critical velocity 
depends on the parameters in the problem. It goes to zero as A ^ and 
the solitons cease to be classically stable. It goes to the speed of light as 
A ^ oo and the barrier to decay becomes infinite. However, Vc does not rise 
particularly rapidly with A: with other parameters chosen as in Fig. ^, Vc 
is only half the speed of light for A ~ 2Ac. Thus, soliton destruction does 
not require that the theory have a value of A lying in some narrow range 
just above Ac- The impact parameter b need not be finely tuned either. Not 
surprisingly, Vc is lowest for collisions with 6 = 0. However, Vc increases by 
less than a factor of two for b of order the soliton radius. Vc also depends on 
the relative orientation angle a between the two solitons in the initial state. 
Here too, the dependence is weak. In the example we explored in detail, we 
found that as a changes from 0° to 180°, Vc varies between 0.25 < Vc < 0.3 
and 0.43 < Vc < 0.48. Thus, although Vc does depend on A and on the 
parameters other than the velocity needed to fully specify a choice of initial 
conditions, the variation of Vc is not dramatic. Soliton decay is not restricted 
to specially chosen velocities, impact parameters, orientations, or values of 
A. Soliton decay is a generic outcome of soliton-soliton collisions. 

Our findings motivate future investigation of collisions between 
metastable solitons in the 3 + 1-dimensional electroweak theory. Previous 
work on two-particle collisions involving these electroweak solitons has 
focussed on collisions between a W boson and a soliton [0. In such collisions, 
the probability for soliton decay falls exponentially as the (rough) analogue 
of A is increased above the (rough) analogue of Ac. This was traced to 
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two facts: First, causing one of these solitons to decay requires delivering 
sufficient energy to one particular mode of oscillation of the soliton. Second, 
a generic incident H^-boson couples very weakly to the mode which must be 
energized if decay is to be induced. We find no analogue of this difficulty 
in our analysis of soliton-soliton collisions in 2 + 1 dimensions. If there 
is a particular mode which must be excited, then soliton-soliton collisions 
seem to generically deliver energy to this mode. And, we certainly see no 
evidence of soliton decay being restricted to theories with |A — Ad <^ Ac. 
This suggests that collisions between two TeV scale particles which can be 
modeled as electroweak solitons (rather than between one VT-boson and one 
such particle) may be an arena in which two-particle collisions generically lead 
to baryon number violation. As we stressed in the Introduction, however, the 
metastable baby skyrmions we analyze differ in several important qualitative 
respects from metastable electroweak solitons. Furthermore, our analysis 
has been purely classical whereas the analysis of VT-soliton collisions in Ref. 
is quantum mechanical. Although our results motivate an analysis of 
collisions between electroweak solitons, they should not be taken to provide 
even qualitative guidance as to the outcome of such a study. 
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